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We present some results concerning the Bose-Einstein effect in asymmetric sources in Monte 
Carlo generators. A comparison of LEP data, standard JETSET predictions and results from 
the weight method is given. Possible generalization of the weight factors to an asymmetric 
form is considered. 
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1 Introductory remarks and data 



Recently we observe a renewal of interest in analysing the space-time structure of sources in 
multiparticle production by means of Bose-Einstein (BE) interference. Such analysis followed 
the example of astrophysical investigations of Hanbury-Brown and Twissu. By improving the 
standard approach! it became possible to model this effect in Monte Carlo generators: as the 
"afterburner" for which the original MC provides a source by shifting the momenta H or by 
adding weights to generated events IB. 

The analysis of BE effect in 3 dimensions is supposed to reflect the spatial source asymmetry. 
Such analysis was done for the LEP data at the Z° peak! which have very high statistics and 
good accuracy. In the following we concentrate our attention on the L3 dataia, as the DELPHI 
data0 are parametrized with only two radii, and the OPAL data0 use the like/unlike ratio 
which requires a cut off of the resonance affected regions even in double ratios. 

As in the L3 paperi we use for each pair of identical pions three components of the invariant 
Q 2 = -(Pi -P2) 2 : QhQlutiQlide defined in the LCMS, where the sum of three - vector mo- 
menta is perpendicular to the thrust axis. Similarly we define a "double ratio" using a reference 
sample from mixed events: 

u I \ Pi 1 Pi 
R2{Pl,P2) = -=3=1 



n mix I mix,MC ' 
Pi P2 

This "double ratio" is parametrized by 

R2(Ql, Qout, Qside) = 7[1 + &Ql + eQout + (Qside] ' 

• [1 + \exp(-R 2 L Q\ - RlutQlut - RtideQlide ~ ^PLfiutRhRoutQhQout)] 

The first bracket reflects possible traces of long-distance correlations; the last term in the second 
bracket seems to be negligible when fitting data and will be omitted in the following. 

By fitting the parameters Rl and R s ide we get some information on the geometric radii in 
the longitudinal and transverse directions (respective to the thrust axis). R ou t reflects both the 
spatial extension and time duration of the emission process. 

In the L3 data the fit region in all three variables extends to 1.04 GeV and is divided into 13 
bins, which gives 2197 points fitted with 8 parameters. The fit parameters 5, e and £ are rather 
small; this means the observed BE enhancement is rather well approximated with a Gaussian. 

The fitted values of radii (in fm) are as follows: 

R L = 0.74 ± 0.02ig;8*, R out = 0.53 ± 0.02±°;°5, R side = 0.59 ± 0.01±^ 

We see clear evidence for source elongation: R s id e /RL is smaller than one by more than four 
standard deviations. 



2 Asymmetric effects from symmetric models 

The geometric interpretation of data requires a comparison with the results from the standard 
MC procedures modelling the BE effect. In the L3 paper such an analysis is given for the 
standard LUBOEI procedure built into the JETSET Monte Carlo generator. This procedure 
modifies the final state by a shift of momenta for each pair of identical pions. The shift is 
calculated to enhance low values of Q 2 and to reproduce the experimental ratio in this variable. 
The superposition of the procedure for all the pairs and subsequent rescaling (restoring the 
energy conservation) makes the connection between the parameters of the shift and the resulting 
ratio rather indirect. 

Using the JETSET parameters adjusted to all the L3 data and the LUBOEI parameters 
fitted to describe the BE ratio in Q 2 the authors of the L3 paper calculated the same quantities 



as measured in the experiment. The projections of R2 look qualitatively very similar to the 
experimental ones. However, the fit to the 3-dimensional distribution gives results different 
from data. The ratio R s id e /RL is not smaller but greater than one; the fitted values (in fm) are: 



R L = 0.71 ± 0.01, R out = 0.58 ± 0.01, R 



side 



0.75 ±0.01. 



We confirmed these numbers in our calculations. We found also that the results are sensitive 
to the JETSET parameters. Using the default values instead of the L3 values we obtained a 
significantly smaller value of R ou t (below 0.5) and significantly smaller A. Other values are less 
affected and R S id e /Ri still bigger than 1. 

We have checked also how the results depend on the source radius R and incoherence pa- 
rameter Xi n assumed in the LUBOEI input function Rbe(Q) = 1 + Xi n exp{— R 2 Q 2 ). In all cases 
we get Rside > Rl > Rout, although the input function was obviously symmetric. The values 
of R s ide and Rl are proportional to R, whereas R ou t changes much less; the dependence on Aj„ 
is very weak. The output value of A decreases quite strongly with increasing R and increases 
with Xi n . No choice of input parameters gives the values of Ri compatible with data. This is 
shown in Fig. [l|. 

Another interesting observation is that to fit the L3 data one needs A = 1.5, which is beyond 
the physically acceptable value of 1. This supports our doubts about the usefulness of the 
LUBOEI procedure in understanding the experimental results (although certainly it is the most 
practical description of data). 
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Figure 1: Fit parameters A and Ri as functions of the input parameters a) for the LUBOEI procedure, b) for the 
weight method. Experimental values are shown on a separate vertical axis. 

Therefore we have also compareclthe data with the results from another procedure modelling 
the BE effect - the weight method 0. In this method each event gets a weight calculated by 
summing the products of two particle weight factors, which are just 1 for equal momenta and 



vanish for large separation in momentum space. A reasonable description of the effect in Q 2 is 
obtained with a simple gaussian form of the weight factor 

w 2 {p, q) = exp[-(p - q) 2 Rl/2], (1) 

or, even simpler, 9 - function form with W2 = 1 for some range of — (p — q) 2 < 1/Rq and ^2 = 
outside. 

In this method we may repeat the same calculation as done for the LUBOEI procedure. 
The resulting double ratios are not that smooth and monotically decreasing as in the data or 
from the LUBOEI procedure (which is the usual drawback of the weight methods). However, 
the major features are surprisingly similar: with weight factors depending only on Q 2 we get 
different values of fitted Ri parameters. Moreover, the hierarchy of parameters is the same: 
Rside > Rl- This suggests that the assymetry is generated by the jet-like structure of final 
states and not by any specific features of the procedure modelling the BE effect. In Fig. lb 
we show the values of the fit parameters as functions of Rq for a Gaussian as well as the 0-likc 
weight factors. Again, no choice of input parameters allows to describe the data. 

These results suggest also that one should be careful with the geometric interpretation of 
the data. If one gets asymmetric distributions from the generator without assuming explicitly 
space asymmetry of the source, it is not clear how the assumed asymmetry will be reflected in 
the results. 

We tried to get some information on this problem within the asymmetric weight method, 
i.e. introducing weight factors which depend in a different way on Ql, Q s id e and Q ou t- This 
work is progress. However, it seems rather difficult to reproduce the data even with two more 
parameters. 



3 Conclusions and outlook 



In this note we investigated the asymmetry of the BE effect in two procedures imitating this effect 
in the Monte Carlo generators and compared them to the data at Z° peak. Both procedures 
give surprisingly similar results and disagree with data. The work on the possible introduction 
of asymmetry in the weight method is in progress. 
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